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Abstract. Resonating valence bond (RVB) liquids in two dimensions are believed 
to exhibit topological order and to admit no local order parameter of any kind. This 
is a defining property of "liquids" but it has been explicitly confirmed only in a few 
exactly solvable models. In this paper, we investigate the quantum dimer model on 
the triangular lattice. It possesses an RVB-type liquid phase, however, for which the 
absence of a local order parameter has not been proved. We examine the question 
numerically with a measure based on reduced density matrices. We find a scaling of 
the measure which strongly supports the absence of any local order parameter. 
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1. Introduction 

The existence of short-ranged resonating valence bond (RVB) liquids has been proposed 
in several two-dimensional quantum spin models [lj. These liquids exhibit no kind 
of simple order but the ground states (GS) are degenerate if the lattice is put on a 
surface with a non-trivial topology (cylinder, torus, etc.). It is (widely) believed that 
such degeneracy has a purely topological origin and is not ascribed to spontaneous 
symmetry breaking in terms of any local order parameter. This topological nature 
has been rigorously confirmed in a few exactly solvable models [2, [3], but proving this 
property in more general cases is a challenging issue. 

In this paper, we address this issue by formulating the problem in terms of reduced 
density matrices (RDM) as we proposed in Ref. [3]. In this formulation, we examine 
the RDMs of the degenerate GSs for various subareas of the system and search for an 
operator distinguishing the degenerate GSs, which means that it can be used as an order 
parameter. We define a convenient measure for this purpose, which is non-zero for an 
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area where an order parameter can be defined. In conventional orders, the measure is 
non-zero on some finite local area, which indicates the existence of an order parameter 
on that area. In RVB liquids, in contrast, it is expected that the measure is zero (in the 
thermodynamic limit) on any local area, thereby verifying topological order. 

With this formulation, we analyze the quantum dimer model (QDM) on the 
triangular lattice, which is one of the simplest microscopic models realizing RVB liquids 
[I]. We numerically show that the dimer liquid in this model cannot be characterized 
by any local order parameter. 

2. Formulation 

2.1. General Case 

We formulate a method to detect an order parameter starting from a general setting. 
Let q be the degeneracy of the GSs and (i — 1, ■ ■ • , q) be the orthonormal GSs. For 
a state \^) in the ground-state subspace, we define its RDM on an area Q by tracing 
out the degrees of freedom outside Q: pq = Tr n where Cl is the complement of 

Q. We also introduce as a reference the RDM averaged over the ground-state subspace: 
^ref = i rjn_ YH=i Note that p T Q is independent of the choice of the basis {|$t)} 

of the ground-state subspace. An order parameter can be defined on Q if and only if 
there exists |^) such that pq ^ p T Q . 

To quantify to what extent RDMs are distinguishable, we introduce a measure Dq 
as follows. We first define a measure of difference between two RDMs as [3] 



where Oq is a (variational) Hermitian operator on Q whose norm is less than unity, i.e., 
|( , 0|C , q|?/>)| < 1 for any normalized vector \ip). Using the eigenvalues {Xj} ofp^— p^ f , this 
can be simplified as: diff (pn,p^ f ) = |A 3 -|. We define Dq by maximizing diff (p^, pff f ) 
over the choice of a GS |^): 



The value of Dq in the thermodynamic limit tells us whether an order parameter can be 
defined on Q. The measure Du has the following useful properties: (a) Normalization to 
a definite range < Dq < 2 — 2/q. (b) Monotonicity: if an area A completely contains 
an area Q, we have Dq < D\. 

How to calculate Dq is not clear from the definition. In the case of QDMs, a 
simple formula is available, which we derive in the next subsection. However, for the 
completeness of the method, here we briefly describe a general algorithm called iterative 
maximization. The difficulty in calculating Dq resides in the double maximizations in 
Eqs. [1] and [21 but we notice that each maximization is possible if either Oq or is 
fixed. The idea is to alternately perform these two kinds of maximizations starting from 
a random vector This would lead at least to a local maximum in the space of Oq 



diff(pn,P?f) = 



max Ttq(Oq P q)-Tiq(Oqp^) , 



(1) 



Dq = max diff (p n , Pq ). 



(2) 
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and |\&). Starting from various random vectors, we would obtain the global maximum 
This method has a potential to handle a wide range of problems. 

2.2. Case of Quantum Dimer Models 

We limit the setting to the case of gapped dimer liquids in QDMs on two-dimensional 
lattices, where further argument is possible. The gapped dimer liquids have been found 
in several QDMs [HE] and related models [21 Ej- However, except for the special cases 
[21 El E] where the correlation length is strictly zero, it is difficult to prove the absence 
of all the local order parameter. This motivates us to study more general cases. 

Let S be the set of all the dimer coverings of the lattice. If the lattice has a non- 
trivial topology, S can be grouped into topological sectors which are not mixed by any 
local operation. We concentrate on the case of the torus in the following. We draw 
two incontractible loops Ai and A 2 which pass through the bonds and wind the torus 
in x and y directions respectively. We classify S into four topological sectors S p with 

p = ++, H — , — h, , depending on the parity of the number of dimers crossing A x 

and A 2 . In a gapped dimer liquid, the lowest-energy states |$ p ) in different topological 
sectors S p become degenerate in the thermodynamic limit. We employ |$ p )'s as the 
basis of the ground-state subspace. 

To define RDMs for QDMs, we must specify the local degrees of freedom of the 
models. To this end, we assign an Ising variable o~k to each bond k of the lattice as in 
Ref. [H] and identify the presence/absence of a dimer on the bond as = +1 and — 1, 
respectively. Any physical configuration {cr^} must satisfy the hard-core constraints: 
for each site of the lattice, there must be exactly one bond with a^ = \ emanating from 
it. An area Q is defined as a set of bonds. We define the matrix element of the RDM of 
a GS |^) as 

(ci\pn\c 2 ) = ^(ci,c|*)(*|ca,c), (3) 

c 

where c± and c 2 are dimer configurations on Q and the sum is over all the dimer 
configurations c on 0. Note that we set (c, c\^f) = if (c, c) is an unphysical configuration 
(violating the hard-core constraint). 

We proceed to the evaluation of D^. We first assume that the area Q is topologically 
non-trivial, i.e., Q encircles the torus. Without loss of generality, we assume that Q 
contains Ai. Let us consider the following trial GS and trial operator: 

l* (trial) ) = !*++>, = J2\c)P(c)(c\, 

c 

where in the second equation c runs over all the dimer configurations on Q and P(c) is 
its parity along A^ These gives a lower bound of D^. 

D n > C>[f ial) (Tr n |*( triaJ )) ( trial ) | - p£) 

Thus the value of is consistent with the existence of a non-local order parameter. 
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We next assume that Q is a (finite) local area. In this case, we can prove the 
following relation: 

Tr n |$ p )($ p ,| = (p^p'). (4) 
To prove this relation, we explicitly express the matrix element of the LHS: 

(d | ( Tin \%) (%>\) |C2> = ($ 

c 

Since f2 is local, we can choose Ai and A 2 so as not to touch Q. Then the dimer 
coverings (ci, c) and (c 2 , c) always belong to a common topological sector and therefore 
(ci,c|$ p ) and (<Ey|c 2 ,c) cannot become non-zero at the same time. Hence we obtain 
Eq. HI Using Eq. HJ we can derive a simpler expression for D^: 

D n = maxdiff(p n ,p£ f ), with p£ = Tr^ |$ p )($ p |. (5) 

p 

This is computable numerically and is employed in the calculation shown below. 
3. Numerical Result 

We consider the QDM on the triangular lattice. The Hamiltonian reads [3]: 

H = E [-*(|//)(jr|+^)+«(|jr)(jr| + |//)(//|)]- (6) 

rhombi 

At the Rokhsar-Kivelson (RK) point (v = t), the GS in each sector is exactly the 
equal-amplitude superposition of all the dimer coverings belonging to that sector [Sill]: 



1 



I*,) = 7pE IO- (7) 

It has been numerically shown that the degeneracy of the GSs and the exponential 
decay of the dimer-dimer correlation at the RK point persist up to some range in the 
parameter space, forming a dimer liquid phase in 0.7 <v/t<l [UEJ Ell ITT] - 

We first numerically calculated -Dn for the RK wave function, Eq. [71 through the 
direct enumeration of dimer coverings (up to N = 52) or the evaluation of Pfaffians 
of Kasteleyn matrices [12] (up to N = 256). The former method is suited for large 
areas while the latter can handle large systems. The lattice is put on the torus and is 
defined by two vectors Ti and T 2 specifying the periodicity. We require that the lattice 
be symmetric under the 27r/3 rotation, choosing Ti — lu + mv, T 2 = —mu + (I + m)v, 
where I and m are integers and u and v are unit vectors as shown in Fig. [TJ The total 
number of sites is given by iV = I 2 + lm + m 2 . As the choice of Q, we define circular areas 
in the following way: we draw a circle centered around a site or an interior of a triangle 
and regard every bond whose midpoint is in the circle as an element of the area (see 
Fig. [T|. The values of Dq, for circular areas are plotted versus the radius R in Fig. [2] 
(a). Though Dq tends to increase as a function of R, it remains small if R is sufficiently 
smaller than the linear system size v~N. It appears that an exponential dependence 
Dq « ce aR ~ b ^ fits the data relatively well, as can be seen in Fig. [2] (b). We also 
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performed calculation for v/t < 1 by obtaining GSs by exact diagonalization. We fit the 
data in the same way and the resultant parameters a, b are shown in the inset of Fig. [3] 
(a). We observe that the relation a ~ 2b holds approximately in the dimer liquid phase. 
This indicates that the finite size effect in Dq enters through y/N — 2R. To interpret 
this, let us consider a loop encircling the torus and passing through fl; see Fig. [3] (b). 
Some part of a loop is not covered by O and the minimal length of the part is VN — 2R. 
The lack of this part makes the distinction of sectors obscure, hence the exponential 
reduction in Dq. Our result explicitly shows the topological nature of the liquid since 
Dq goes to zero exponentially with iV — > oo for any fixed R. This does not only exclude 
local order parameters which are diagonal in the dimer basis (and amenable to exact 
calculations in the iV — > oo limit [llj), but any, possibly non-diagonal, operator. The 
parameter b is naturally interpreted as the inverse of a correlation length, and the value 
we obtained at the RK point is close to ~ 0.83 ±0.12z7r and ^ 1 ~ 0.76 obtained by 
Ref . [11] , where A refers to the direction along a bond and B the direction perpendicular 
to it. 
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Figure 2. The result for the RK wave function, (a) The value of Dq as a function of 
the radius R for different system sizes N. (b) Fitting of the data using an exponential 
function Do « ce aR - b ' /W . 
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Figure 3. (a) Main panel: fitting of the data of Dn for v/t = 0.9. Inset: obtained 
fitting parameters a, b for different values of v/t in the dimer liquid phase, (b) A loop 
(broken line) encircling the torus and the area fi. Periodic boundary conditions are 
imposed along T x and T 2 . A part in the loop with the length of \/N — 2R is not covered 
by 0. 



4. Conclusion 

We formulated a method to detect the existence of an order parameter using reduced 
density matrices. We numerically applied the method to the liquid phase of the QDM 
on the triangular lattice. The measure Dq we defined contrasts local areas with non- 
local areas: Dq — > with N — > oo on local areas while Dq > 1 on non-local areas. The 
data of Dq is fitted well by an exponential function, which can be interpreted from the 
topological picture. Our result explicitly verifies topological order in this system. 
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